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ABSTRACT 

While there are numerous indications that GRBs arise from the death of massive stars, the GRB 
rate does not follow the global cosmic star formation rate and, within their hosts, GRBs are more 
concentrated in regions of very high star formation. We explain both puzzles here. Using the publicly 
available VESPA database of SDSS Data Release 7 spectra, we explore a multi-parameter space in 
galaxy properties, like stellar mass, metallicity, dust etc. to find the sub-set of galaxies that reproduce 
the recently obtained GRB rate by Wanderman & Piran (2010). We find that only galaxies with 
present stellar masses below < 10 10 M© and low metallicity reproduce the observed GRB rate. This 
is consistent with direct observations of GRB hosts and provides an independent confirmation of the 
nature of GRB hosts. Because of the significantly larger sample of SDSS galaxies, we compute their 
correlation function and show that they are anti-biased with respect to the dark matter: they are in 
filaments and voids. Using recent observations of massive stars in local dwarfs we show how the fact 
that GRB hosts galaxies are dwarfs can explain the observation that GRBs are more concentrated 
in regions of high star formation than SNe. Finally we explain these results using new theoretical 
advances in the field of star formation. 

Subject headings: galaxies: evolution, galaxies: statistics, galaxies: stellar content, gamma-rays: bursts 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are short and intense 
pulses of soft gamma-rays. It is generally accepted that 
long-duration (> 2 s, but see Bromberg et al. 2012) arise 
durin g the Collapse of massive star s the so called Collap- 
sars (MacFadyen & Wooslcy 1999). This understanding 
is largely based on the observed association of several 
long bursts w ith type Ibc supernovae (SNe) (Hjorth & 
Bloom||2012| Fl In such a case, it is only natural to ex- 
pect that GRBs will follow the cosmic star formation 
rate (SFR) and that like SNe, their positions within their 
host galaxies will correspond to the star formation rate 
within these galaxies. Surprisingly both expectations are 
not satisfied by observations. 

During the last eight years the Swift satellite has rou- 
tinely provided us with accurately localised GRBs. From 
this data it is possible to construct t he luminosity func- 
tion and c osmic GRB rate. Recently |Wanderman fc PT 



ran 



( 2010[ ) (WP hereafter) have estimated the rate and 
luminosity function of long duration GRBs using a novel 
method that solves simultaneously for the GRB rate and 
the GRB luminosity function. One of their most interest- 
ing findings was that, assuming that the GRB luminosity 
function does not depend on cosmic time, the GRB event 
rate does not follow the star- format ion rate of the typical 
galaxy population (their Fig. 9) , showing deviations both 
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at low (z < 3) and large (z > 3) redshifts. On the other 
hand, using HST ima ging of GRB and SN host galaxies 
Fruchter et all ( 2006 ) showed that while the likelihood 
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to find a SJNe m a given position within its host galaxy 
is linearly proportional to the SFR at that point, GRBs 
are much more conce ntrated within the r egions of the 
highest SFR (see also |Svensson et a"T]|2010 |. 

An obvious question that arises is how these two is- 
sues can be reconciled with the idea that long GRBs 
arise from the collapse of a massive star. A related ques- 
tion is of course what kind of galaxies hosts GRBs and 
whether the SFR within this hosts is similar to the cos- 
mic one. Ample information is available on the nature 
of GRB hosts. Inspection of the hosts ( Fruchter et al 
[19991 |Chary et al.|[2002| |Bloom et al.||2002||Le l'loc4i et 



al.||20031 ITanvir et all|2004| |Fruchter et al.[|2006| |Cas- 
' •• 2006[TSavaglio et al.||2009| l~(see also 



tro Ceron et al . 
Fynbo et al.|2TjT2" 



Levcsque 2013 tor recent reviews) re 



veals that usually they a re low mass irregular g alaxies 
and have low metallicity ( Prochaska et al||2004[ [S oller- 
man et al.|2"005l|Fruchter et al.|2006||Modjaz et al.|2006[ 



Stanek et al||2006| |Thone et al.||20071 |Wiersema et al 



20071 [ Margutt i et al.||2007| |Savaglio et al.||2009[ JThone 
et al. ||20131. These findings are consistent with theoreti- 



cal modeling that suggests that low metallicity is essen- 
tial to produce high angular mome ntum and high-stellar 
mass needed for GRB pr o genitors flYoon fc Langer | 2005 [ 

Podsiadlowskil |2007J). 
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Woosley & Heger 

Note, however, that recent observations have revealed 



a population of more massive G RB hosts (Castro Ceron 
et al|[2006 Savaglio et al. 2009), in particular w hen the 



dustiest galaxies are targeted (Perley et al. 20131. Never 
the less, the fraction of such dusty and massive galaxies is 
not high enough to reconcile the GRB rate with the typi- 
cal SFR, at least out to z=l ( Perley et al|2013 Kocevski 
|et al.|2009| , i.e. GRBs do prefer galaxies with lower mass 
with respect to the typical star-forming galaxy popula- 
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tion. As for the metallicity, while GRB hosts typically 
have low metallicity, a h andful of GRB absorbers show 
roughly solar abundanc e (Savaglio et al. 12003 
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|De (Jia et al.||2012| ). 
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In an attempt to address the first question we ignore 
all this information and attack this question using a dif- 
ferent approach. Our way to answer the question - which 
galaxies host GRB - is to have a complete census of star 
formation histories of all galaxies in the universe and ex- 
tract the ones that match the observed GRB rate. Here 
we do it so by exploiting a large SDSS local sample of 
galaxies and analysing in detail their fossil record. This 
allows us to effectively have a complete sample of galax- 
ies up to z ~ 3 and compare the star formation rate in 
low mass low metallicity galaxies with the inferred GRB 
rate. 

Once we address this question, we turn to the second 
question: why are GRBs more concentrated in high SFR 
rate regions than SNe? To address this issue we make 
the reasonable assumpti on that GRBs arise from more 
massive stars than SNe (Ostlin et al. 2008). Following 



this assumption we compare the positions, within low 
mass galaxies, of massive (> 20Mq) stars with the posi- 
tions of less massive (> 9M Q ) stars. Finally, we suggest 
a simple theoretical argument concerning star formation 
that explains this trend. 

2. METHODOLOGY 

The spectra of galaxies encodes information about the 
histories of the component stellar populations, dust, and 
star formation. Various tools have been deyeloped to 
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extract this inf ormation (e.g., Heavens et al i 
jeiro et aL||2007 ) and previous works have compared the 
resulting extracted information with bot h extrinsic and 
intrin sic galaxy properties. The MOPED (Heavens et al. 
20001 algorithm implements the general process of re- 
forming a complex dataset (e.g., a galaxy spectra) into 
a set of parameters (e.g., star formation rate, metallic- 
ity) and parameter combinations, assuming uncorrelated 
noise, such that the data compression is loss less. 

An easily accessible, robust code, is the VErsatile 
SPectral An alysia^] (hereafter VESPA, see To jeiro et al. 
2007 2009| for more details) package, which, recovers 
star formation and metallicity histories of the galactic 
spectra using synthetic stellar-population models. The 
software recovers histories in adaptive age bins accord- 
ing to the signal-to-noise of the galaxy spectrum on a 
case by case basis and addresses the age-metallicity re- 
lation. Two popular synthetic stellar population models 
are included in the VES PA output ; those of IBruzual ~El 
Chariot! f2003"| , and the |Maraston| ( |2005[ ) fc |Maraston 
et al.| ( |2009| ), which differ in their respective resolutions, 
and the use of empirical libraries to model the thermally 
pulsating asymptotic giant branch. Furthermore VESPA 
corrects for Galactic extinction using the dust maps of 
Schlegel et al. ( 1998 ), and fits for the dust in each galaxy 
using either a one or two parameter dust model. We 
exploit the VESPA database to compare the observed 
GRB rate with that of SDSS galaxies. In particular, we 
address the question: is there any subset of galaxies from 

5 http:/ /www-wfau. roe. ac.uk/vespa/ 
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Table 1 

The 16 time bins within which VESPA determines the star 
formation fraction in the rest frame of the galaxy, which 
corresponds to the time t = Tb ■ 



the SDSS that matches the observed GRB rate? If so, 
what are the physical properties of these galaxies? 

3. RESULTS 

The ~ 10 6 spectroscopically selected galaxies in this 
work were drawn from the Sloan Digital Sky Survey Data 
Release 7 ( |York et al.1|2000| |Abazajian et al.|l2009[ and 
references therein, hereafter SDSS DR7). We selected 
both "Main Galaxy Sample" and Luminous Red Galaxies 
(LRG). Galaxy s pectra were re p rocess ed using VESPA 
and wc followed 



Tojeiro et al. (2009| in adopting the 



two parameter dust model, but note that our results are 
insensitive to the choice of dust models. 

We are concerned with the star formation histories 
of galaxies. We extract this quantity from the VESPA 
database for the highest resolution output (16 bins) in 
lookback-time bins as specified in Table 1. VESPA also 
returns the total stellar mass, both present and at for- 
mation, of the galaxy and the metallicity for each bin. 
We explore total mass and metallicity in the VESPA 
database in order to see if a match is found to the WP 
data (their fig. 2, upper panel). Note that the derived 
WP GRB rate shows two characteristic features that do 
not match the global star formation history of the gen- 
eral population: the peak is shifted to z ~ 3 and below 
that redshift the GRB rate distribution is flatter than the 
SFR. At high red shift the GRB rate is higher than the 
SFR rate (see also |Jakobsson et al. 2012 ). From previous 
experience, we know that this beha viour is a character- 
istic of low-mas s gala xies ( Heavens et al.|2 004). Fig. 2 in 
Heavens et al.| |2004[ ) shows how the shape and peak of 
the star formation history of galaxy changes as a func- 
tion of mass; this already points out to low-mass galaxies 
as promising candidates. This idea is also supported by 
the non-detection of extreme-redshift GRB host galaxies 
in very deep HST imaging (Tanvir et al. 2012). How- 
ever, in order to find a good match to the WP GRB rate 
a cut in metallicity is also required. In particular, we 
find that if we select SDSS galaxies with present stellar 
masses < 10 10 M Q and metallicity Z < 0.1Z Q , we obtain 
a good match between the GRB rate and the SFR. This 
means that we need to exclude about half of the dwarfs 



to obtain such metallicity cut-off (see Fig. 6 in Panter et 



log(M/M s ) 



Figure 1. The co-moving rate of GRBs from WP (diamonds) with 
the best fitting star formation history obtained from exploration of 
the VESPA catalog. Only after imposing a cut-off in galaxy stellar 
mass (< 10 10 Mq) and metallicity (< 1/10Zq) one finds a good 
fit (solid line). 



al.||2008 ) . Any other population produces a steeper SFR 
with redshift and thus cannot be matched to the GRB 
rate. 

Our reported metallicities correspond to the ones of 
the stellar population in the whole galaxy, as derived 
by VESPA. These are different from those typical ly in- 
ferred locally for GRB hosts. For example, recently, |Gra-| 
ham fc Fruchter| ( |2012[ ) have measured the metallically of 
the star forming regions (HII regions) in identified GRB 
hosts. They measure oxygen abundances from emission 
lines. Their derived metallicities are somewhat larger 
than ours (by a factor of 50%). This is due to the fact 
that we measure the whole stellar population, including 
also the old stellar population in the irregulars, this low- 
ers the overall me tallicity (see the blue line of Fig. 3 in 
Panter et al.|2008 ). So reporting only the metallically of 
the latest and most recent burst in star formation would 
bring our metallicity in agreement wit h the metallicity 
reported by Graham & Fruchter (2012 1. 

Fig. 1 depicts the WP GRB rate distribution with 
redshift, where the solid line marks the result of of our 
fit. First, note that we could not go beyond z > 3 as the 
VESPA catalog does not contain information on these 
high-redshift bins. The reason for this is that we are 
reconstructing the star formation history from the fossil 
record at a medium redshift of z ~ 0.2 and going back 
so far in time is very difficult. Therefore we limit our 
analysis to redshifts below 3, not covering the peak of 
the observed GRB rate. We also note that at low— z the 
VESPA reconstructed SFR is slightly higher that a the 
WP data at z = 0.25, showing a recent bump of star 
formation activity. 

We find that, based just on matching the observed SFR 
in the specific population of galaxies with the WP GRB 
rate, the hosts of long GRBs are dwarf galaxies with very 
low metallicities. This result is in good agreement with 
direct observations of those hosts ( |Savaglio et al. 2009), 
but these samples are limited to a tew objects (just 46) . 
In fact, the inconsistency between the GRB rate and typ- 
ical star formation could have been predicted, given that 
most GRB hosts are low-metallicity dwarf galaxies with 
a different SFR than the typical star-forming galaxies 



Figure 2. Number of stars as a function of the dwarf galaxy 
mass for five dwarf galaxies. The diamonds correspond to stars 
with masses > 20 Mq while triangles to > 9 Mq. The lines are 
fits to the data points and show that the more massive dwarfs have 
more massive stars (> 20 Mq) with respect to less massive stars 
(> 9 Mq) than lower mass dwarfs. 



and hence the overall SFR. The agreement between the 
metallicity cut in our analysis and the observational evi- 
dence that most GRB hosts have low metallicity further 
supports the WP results on the GRB rate and suggests 
that the specific nature of GRB hosts rather than a lumi- 
nosity evolution, for example, is responsible for observed 
distribution of GRBs redshifts and peak fluxes. 

We now turn to the second puzzle, concerning the ori- 
gin of the discrepancy between the GRB and SN posi- 
tions within their hosts. To do so we explore the proper- 
ties of massive stars in these dwarf low-metallicity galax- 
ies. Specifically, we compare the population of massive 
stars > 9 M with the one of very massive stars > 20 
Mq. The basic idea is that while the former lead to SNe, 
the latter are GRB progenitors and if the two populations 
are distributed differently this explains the different dis- 
tribution of GRBs vs. SNe. 

The most detailed study of the population of massive 
stars i n dwa rfs in the local group is that by |Bianchi| 
ft <i.L| j2012l. who used (tALEX obscry<it ions to deter- 



mine the number of massive stars (> 9 M ) in six dwarf 
galaxies (Phoenix, Pegasus, Sextans A and B, WLM and 
NGC6822), although for the last one they only have lower 
limits. We have used the masses for these dwarfs de- 



rived by Karachentsev et al. (20041 (column 4 in their 



Table 4) to plot in Fig the numbers of massive stars as 



a function of the dwarfgalaxy mass. Inspection of Fig. [2] 
reveals that there is a clear correlation between the num- 
ber of massive stars and the mass of the dwarf galaxy. 
Additionally, the ratio (# > 20)/(# > 9) of number of 
very massive to massive stars increases with the mass of 
th e dwarf galaxy. The li nes show fits to the data points 



Bianchi et al. ( 2012 ) and are well fitted by power 
laws with power 0.7 (for > 9 M ) and 1.2 (for > 20 
Mq), i.e. there is a nearly factor two difference in the 
power. This empirical observation implies that very mas- 
sive stars (> 20 M ) are more abundant in more massive 
dwarfs. This plot can b e compared directly with Fig. 4 



Svensson et al. 



but for number o 



(2010), which plots the same quantity 
SJN and GRBs in hosts galaxie s . Our 
results agree very well with the |Svensson et aL (20101 




current consensus (McKee & Ostriker|2007 ) is that mas 



Figure 3. Number of stars with mass > 9 Mq (dashed line) or 
> 20 Mq (dotted line) as a function of the total light in the dw arf 
galaxy compared with the results from iFruchter et al.l (120061 for 
the numbers of SNe (upper solid line) and Ciittis (lower solid line), 
note the good agreement. 



trend. Given that larger dwarf galaxies are denser than 
lower mass ones, we can take this relation as a proxy to 
the relation according to which very massive stars are 
more abundant in denser regions in such galaxies. This 
assertion, that we will shortly explore further, serves now 
as the basis for the rest of the analysis. 

We can use the data in |Bianchi et al.| (|2 012|_ to pro - 
duce a plot similar to Fig. 3 in |Fruchter et al.| (J2006), 
which shows the number of SN/GRBs as a function of 
the total light in the galaxy. Because we do not have 
pixel information for the dwarfs, we simply integrate the 
light in radial r i ngs. T his is equivalent to what is done n 
Fruchter et al.| ( |2006[ ) if the galaxies do not show angu- 
lar asymmetries, which is the case for the dwarfs we are 
considering. The results are shown in Fi g . [3} T he solid 
thick lines are fits to the Fruchter et al. (2006) results, 
while the dashed line is the number of stars > 9 M Q and 
the dotted line > 20 Mq, where we have averaged over 
the five local dwarfs: there is good agreement. 

Because dwarf galaxies have very low surface density, 
the outer edges of dwarf galaxies will co ntain very lit- 
tle m ass as compared with the inner part (Tolstoy et al. 
2009 1 . For concreteness let us denote as outer edges' 
all mass beyond one equivalent radius and as the "inner 
part' all mass inside it. For dwarf galaxies this means 
that there is less than a factor 10 mass in the outer edges 
than in the inner part. If we now use Fig. [2] this means 
that the outer edges will have a factor 3 less very massive 
(> 20 M Q ) stars relative to massive stars (> 9M ) than 
the inner part. If we now assume that GRBs can only 
originate on very massive stars then this is in agreement 
with the obse rvations of the spatia l distribution of SNe 
and GRBs by |Fruchter et al.| (|2006| ). Thus the fact that 
GRBs are mostly hosted in dwarf galaxies, when com- 
bined with the dependence of the initial mass function 
on the density in these galaxies, explains also the spatial 
distribution discrepancy between SNe that follow the lo- 
cal SFR linearly and GRBs that follow the local SFR 
much stronger than linearly. 

It is interesting to explore whether there are any theo- 
retical hints that would suggest this observed "Fruchter 
et al. law" of a differential initial mass function. The 



sive stars form beca use of the turbulent nature of the 
interstellar medium. Padoan & Nordlund (2002) give a 
detail theory of star formation due to magnetised turbu- 
lence. In particular they predict that the maximal mass 
of the initial mass function is (their eq. 16): 



pL 3 



(1) 



where p is the density in the star forming region, L the 
size of the largest scale where turbulence is driven and 
A4a the Alfven Mach number. It is clear from the above 
equation that lower density, implies a lower maximum 
mass. This is the case for dwarf galaxies, in which the 
density beyond the equival ent radius, is a fact or few be- 
low that in the inner part ( Tolstoy et al.|2009 ). The size 
of the largest turbulent scale will remain tfie same while 
the M.a will be similar or slightly larger because of the 
lower density. Thus, the maximal mass in the outer lower 
mass region will be lower as observed (see Fig. [2]). Note 
also that the same authors predict that at lower (Tensities 
the peak of the mass function will move to larger masses, 
thus increasing the number of SN. 

While the above theoretical argument and the observa- 
tions of massive stars in local dwarf galaxies would sug- 
gest that the number of GRBs should equal the number 
of SNe as the mass of the host galaxy increases beyond 
10 10 Mq, it is the metallicity cut-of that prevents this 
from happening. 

4. IMPLICATIONS 
The first implication of our results is that, as su ggested 



by direct observations ( Graham fc Fruchter 1 20 121) and by 
theoretical arguments ( Yoon & Langcr 200o| |\Voosley &; 



Heger|2006||Wolf fc Pod siadlowski 2007|), low metallically 
is an important part of the Collapsar model. While some 
special cases show GRBs in high metallicity regions, low 
metallicity is clearly an important factor for most of the 
population. This has numerous implications on models 
of GRBs' progenitors and on the operation of their inner 
engines. 

The fact that low-metallicity plays an important role in 
the formation of GRBs has also far reaching implication 
concerning the use of GRBs to explore the early universe. 
Properly utilising GRBs as probes of the early universe 
requires a thorough understanding of thei r formation and 
the host environments that they sample ( |Levesque|2013 ). 
Galaxies exhi bit a strong mass - metallicity relation (see 
e.g. Fig. 6 in Panter et al.|2008 ) and massive galaxies (> 
10 11 Mq) do not have such low metallicities as those of 
GRB hosts. Thus, the host galaxy population of GRBs is 
biased relative to the overall galaxy population and even 
relative to the population of dwarf galaxies as we had to 
impose a metallicity cut of w 1/10 the solar value, this 
will exclude about half of the dwarf galaxies as possible 
hosts of GRBs. 

This last point suggests that as we move to the early 
universe, where metallicity was lower, the GRB rate was 
much higher and a significantly larger fraction of stars 
resulted in GRBs. This is consistent with the observa- 
tion in WP that the GRB rate at high redshifts is signifi- 
cantly flatter than the SFR. Note however, that a massive 
galaxy will enrich its gas extremely fast (in a dynamical 
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Figure 4. Sky location for the host galaxies of visible GRB, that are needed to match the WP GRB rate, is denoted as blue dots. The 
lower slice corresponds to z = 0.1 and the upper slice to z = 3 (assuming the current location). For reference, we also show galaxies 
classified as luminous red galaxies, i.e. massive (> L„) galaxies with mostly old stellar populations. It is apparent that the hosts of visible 
GRBs occupy mostly the filaments and voids of the cosmic web. 

a high metallicity of the available gas to form stars; this 
explains why the GRB rate decreases at lower redshift. 
However, because dwarf galaxies dominate the SFR at 
low-z, the GRB rate will decrease slower than the SFR 
of the overall galaxy population. If we use the metallic- 
ity of the damped Lyman— a systems as a tracer for the 
metallicity of the gas where galaxies form, we see that 
at redshift z w 2 — 3 the metallicity of th e gas decreases 
below 1/10 the solar value (see Fig. 12 in |Rafelski et al. 
( 2012 1). This same conclusion is obtained it one looks at 
the metallicity hist ories of SDSS galaxies (see Fig. 3 in 




r/Mpc 

Figure 5. Correlation function for the star formation rate in the 
hosts of visible GRBs as a function of scale. The grey band shows 
the computed range for the dwarfs in the SDSS DR7 needed to 
reproduce the GRB rate in WP. Note that GRB hosts are less con- 
centrated than the mean population (dotted line), thus indicating 
that their pair separation is larger that the mean of the galaxy 
population. 

time), thus even if the original gas is of low metallicity, 
only galaxies that have low rates of star formation will 
be able to host GRBs as they will be able to keep the gas 
at low metallicity. As we move to lower redshift and the 
overall metallicity of the gas enriches, even dwarf galaxies 
will have difficulties at hosting GRBs if they already have 



Panter et al. 20081, thus the GRB rate should increase 
at that redshift in agreement with WP. 

After having identified GRB hosts as low-metallicity 
dwarf galaxies, we can now explore properties of this 
population within the SDSS data. First, we examine the 
location in the sky as a function of RA and declination of 
GRB hosts. Fig. 14] depicts the GRB hosts as blue dots. 
For reference we have also plotted the locus of luminous 
red galaxies (i.e. galaxies with larger masses (~ L„) with 
older stars in them). It is visually apparent that the hosts 
of GRBs are situated in the filaments- voids of the cosmic 
web. The bottom panel of Fig. [4] shows the positions at 
z = 0.1 while the top panel show their position at z = 3 
(for this case we have used the current local positions but 
identified the galaxies according to their star formation 
history at z = 3). 

In order to quantify the location in the sky of GRB 
hosts we compute their correlation function, weighting 



G 



by the star formation history, thus we compute the so- 
called mark correlation, using star formation as a mark. 
As can be seen in Fig. [5] where the notation WW/DD 
indicates that the mark statistic is the ratio of weighted 
pair counts to unweighted pair counts (in this notation, 
the traditional unweighted correlation function would be 
DD/RR, where RR is the number of un-weighted pair 
counts in a random distribution.). The hosts of GRBs are 
less correlated than the mean population at scales smaller 
than 20 Mpc. The shaded region shows the 2a confidence 
region from all GRB hosts in the SDSS. Our prediction 
then is that GRBs are to be found in under-dense re- 
gions of the dark matter density field, i.e. they will be 
anti-biased. Most GRBs inhabit regions that show the 
lowest rates of merging and are undisturbed, despite ev- 
idenc e that some systems do show signs of interaction 
(e.g.,|Fruchter et al.||2006|). 



Levan, A. J., Strolger, L., et al. 2006, Nature, 



It is curious to point out that planet formation requires 
exactly the opposite environment: high-mctallicity. 
Therefore, GRBs overall seem to inhabit regions where 
no planets form, thus no presenting a risk for life exter- 
mination. While this is generically correct all over the 
universe, it is particularly relevant concerning life exter- 
minations on Earth. The Milky Way is a rather massive 
galaxy with very few low metallicity regions. In fact, us- 
ing the star format ion and metallicity histories of Milky 
Way type galaxies ( Panter et al.|2008| in the SDSS sam- 
ple we find that only 2% of the Milky Way has metallicity 
below of 1/10 of the solar value which we found as typical 
upper limit to galaxy that host GRBs. This implies that 
the expected rate of GRBs in the Milky Way is much 
below the one expected from a simple estimate based on 
the local cosmic GRB rate per unit time and unit volume 
(~ l.SGpc^yr- 1 ). 
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